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ABSTRACT 

Soft extended emission (EE) following initial hard spikes up to 100 seconds was observed with 
Swift/BAT for about half of short-type gamma-ray bursts (SGRBs). This challenges the conversional 
central engine models of SGRBs, i.e., compact star merger models. In the framework of the black 
hole-neutron star merger models, we study the roles of the radial angular momentum transfer in the 
disk and the magnetic barrier around the black hole for the activity of SGRB central engines. We 
show that the radial angular momentum transfer may significantly prolong the lifetime of the accretion 
process and multiple episodes may be switched by the magnetic barrier. Our numerical calculations 
based on the models of the neutrino-dominated accretion flows suggest that the disk mass is critical 
for producing the observed EE. In case of the mass being ~ O.8M , our model can reproduce the 
observed timescale and luminosity of both the main and EE episodes in a reasonable parameter set. 
The predicted luminosity of the EE component is lower than the observed EE with about one order of 
magnitude and the timescale is shorter than 20 seconds if the disk mass being ~ 0.2A/©. Swift/BAT- 
like instruments may be not sensitive enough to detect the EE component in this case. We argue that 
the EE component would be a probe for merger process and disk formation for compact star mergers. 

Subject headings: Gamma-ray burst: general-accretion, accretion disks - black hole physics 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) are sorted into two classes 
(Kouveliotou et al. 1993), i.e., short duration (Tgo < 2s, 
SGRBs) and long duration GRBs (T 90 > 2s). Their pro- 
genitors are thought to be mergers of two compact stars 
(Eichler et al. 1989; Paczyhski 1991; Narayan et al. 1992; 
Zhang et al. 2007; Nakar 2007) and collapses of massive 
stars (Woosley 1993; Paczyhski 1998; Piran 2004; Zhang 
& Meszaros 2004; Woosley & Bloom 2006), respectively. 
However, the observed burst duration is instrumentally 
dependent (Donaghy et al. 2006; Qin et al. 2012). Swift 
observations reveal that the short-long GRB classifica- 
tion scheme does not always match their physical origin 
classification scheme, i.e., mergers of compact binaries 
(Type I) vs. collapsars (Type II) (Zhang 2006; Zhang 
et al. 2007; Zhang et al. 2009; Lu et al. 2010; Xin et 
al. 2011). With the CGR O/BATSE data, Lazzati et al. 
(2001) found an excess emission peaking ~ 30 s after the 
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prompt one, which is detectable for ~ 100 s for some 
SGRBs (see also Connaughton 2002; Norris et al. 2010). 
About half of the lightcurves of those GRBs that are 
recognized as Type I GRBs with Swift/BAT show ini- 
tial hard spikes following by an extended emission (EE) 
component of soft gamma-rays up to ~ 100 s post the 
BAT trigger (Barthelmy et al. 2005; Lin et al. 2008; 
Perley et al. 2009; Zhang et al. 2009). The most promi- 
nent case is GRB 060614. Its lightcurve is composed 
of some initial hard spikes and a long, soft gamma-ray 
tail, which leads to its T90 being ~ 110 s (Gehrels et 
al. 2006). The initial hard spikes could be recognized 
as a SGRB with CORO/BATSE-like instruments, since 
the soft EE is out of the instrument bands (Zhang et 
al. 2007). No accompanied supernovae was detected for 
this nearby long GRB (the redshift z = 0.1254; Delia 
Valle et al. 2006; Fynbo et al. 2006; Gal- Yam et al. 
2006), disfavoring the collapse of a massive star as the 
progenitor of this GRB. On the other hand, it is shown 
that some intrinsically short GRBs are likely of Type 
II origin (Zhang et al. 2009; Belczynski et al. 2010; 
Levesque et al. 2010; Lin et al. 2010; Xin et al. 2011; 
Virgili et al. 2011). These observations indicate that Tgo 
may be not a good parameter to distinguish two types of 
GRBs. A detailed analysis on the instrumental selection 
effect and energy dependence of T90 with Fermi/GBM 
data by Qin et al. (2012) further supports this idea. Lii 
et al. (2010) proposed a new classification parameter, 
i.e., e = Eisa/Epl, to group an observed GRB into the 

physically- motivated Type I/II classification schem^El 
where Ei so is the isotropic gamma-ray energy and E P:Z is 
the peak energy of the vf v spectrum in the rest frame. 
They showed that some SGRBs are sorted into the high-e 
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group as typical Type II GRBs. 

The observed EE component challenges not only the 
short-long GRB classification scheme, but also the con- 
ventional central engine models for SGRBs. Popu- 
lar central engine models of Type I GRBs are related 
to the accretion on to a central compact objects that 
is formed from merger of a stellar compact binary, 
namely neutrino-dominated accretion flows (NDAFs, 
e.g., Popham et al. 1999; Narayan et al. 2001; Di Mat- 
teo et al. 2002; Kohri & Mineshige 2002; Kohri et al. 
2005; Gu et al. 2006; Chen & Beloborodov 2007; Liu 
et al. 2007, 2008, 2010a, 2010b, 2012; Lei et al. 2009; 
Sun et al. 2012). Merger of a black hole (BH)-neutron 
star (NS) binary is the most favorite scenario. Such a 
system would result in a rotating BH with several so- 
lar masses surrounding by a neutrino-cooled disk. The 
detection of the EE component likely suggests that the 
central engine is not died out rapidly. Several lines of 
evidence from Swift/BAT observations also support this 
idea. It was also proposed that the early shallow decay 
X-ray emission and internal X-ray plateau may be due 
to the spin-down energy release of the proto-magnetar of 
a stellar compact binary merger (Dai & Lu 1998; Zhang 
& Meszaros 2002; Lyons et al. 2010). Late X-ray flares 
may signal the restart of the GRB central engine and 
evolution of the disk (Fan & Wei 2005; King et al. 2005; 
Burrows et al. 2006; Dai et al. 2006; Perna et al. 2006; 
Proga & Zhang 2006; Lazzati et al. 2008; Lee et al. 
2009; Yuan & Zhang 2012). Lazzati et al. (2008) inves- 
tigated the temporal evolution of the disk in GRB central 
engines to explain the observed decline of X-ray flare lu- 
minosity. They argued that it is the dynamics of the disk 
or the jet launching mechanism to generate an intrinsi- 
cally unsteady outflow on time-scales much longer than 
the dynamical timescale of the system for the late X-ray 
flares. It was also suggested that propagation instabili- 
ties, rather than variability in the engine luminosity, are 
responsible for some X-ray flares (Lazzati et al. 2011). 

Different from the early shallow decay X-ray emission 
and internal X-ray plateau, the EE component is usu- 
ally highly variable and is usually not clearly separated 
from the burst itself in their lightcurves. It may be pro- 
duced by the same process as the prompt emission in 
different episodes. Mctzger et al. (2008) presented time- 
dependent models of the remnant accretion disks created 
during compact object mergers. They calculated the dy- 
namics near the outer edge of the disk to study the evo- 
lution of the accretion rate at a long timescale (100 s or 
longer). They showed that the late-time accretion can in 
principle provide sufficient energy to power the late time 
activity observed by Swift/BAT from some SGRBs. In 
their models, the majority of the disk mass is in the outer 
edge and the disk becomes advective at the late time. In 
this paper, we focus on the radial angular momentum 
transfer in the disk and the magnetic barrier around the 
BH that may affect the activity of SGRB central engines. 
Based on the NDAF models, we present detailed calcula- 
tions and apply our model to some typical SGRBs with 
detection of the EE component. We describe our model 
in Section 2. Numerical results are shown in Section 3. 
Conclusions and discussion are presented in Section 4. 

2. MODEL 



Assuming that the progenitor of SGRBs is a BH-NS 
binary, the merger of this system would result in a ro- 
tating BH surrounding by a neutrino-cooled disk. We 
focus on the roles of the radial angular momentum trans- 
fer in the disk and the magnetic barrier for prolonging 
the life time of the SGRB central engine (e.g., Proga & 
Zhang 2006). A fraction of the disk matter may carry 
part of the angular momentum of the accretion matter 
to form a radial outflow. The competition of the radial 
angular momentum transfer to the gravity of the central 
BH would increase the timescale of the accretion process 
(see, e.g., Lee & Ramirez- Ruiz 2007; Rosswog 2007; Mct- 
zger et al. 2010). On the other hand, the magnetic field 
accumulated near the horizon of the BH may be strong 
enough to prevent the gravity and accretion process. The 
magnetic barrier would dissipate quickly as the accretion 
rate drops. These processes likely lead to multiple well- 
connected accretion episodes (e.g. Narayan et al. 2003; 
Cao 2011). We illustrate the processes of our model and 
corresponding cartoon lightcurve in Figure 1. Our model 
is elaborated below. 

2.1. Outward angular momentum transfer 

Without considering the mass (energy) and angular 
momentum lost in jet production, the conservations of 
the mass (energy) and angular momentum read (e.g., 
Bardeen 1970; Thorne 1974; Wang et al. 2002), 

M n+ i - M n = (M* - M* +1 )e in , n = M n T n e in , n , (1) 

Jn+1 — Jn = J n ~ = M n T n l[ nyn , (2) 

where M n ( J n ) and M* (J*) are masses (angular momen- 
tums) of the BH and disk, respectively, M n is the mass 
accretion rate, T n is the accretion timescale, and ej nin and 
^in,n are the specific energy and angular momentum at 
inner boundary orbit in the nth episode (n = 1, 2, 3, ...). 
The mass of the disk can be written as 

rr ou t,n 

M* = 2tt / S n rdr, (3) 

J ^in,n 

where r m ^ n and r outi n are the inner and outer boundaries 
of the disk in the nth episode, respectively. The angular 
momentums of the BH and disk can be calculated with 

Jn = (4) 

c 

/r ou t,n 
S n /„rrfr, (5) 

where l n is the specific angular momentum per unit mass. 
It is known that a* n cannot exceed unity (e.g., Janiuk et 
al. 2008). Since supplies all the accretion process, 
we have 

M*=^M n T n , (6) 

n 

where T n is the timescale of the nth epoch, which de- 
pends on the competition between the magnetic flux 
pressure and the gravity to the accreting mass. 
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2.2. Switch of the magnetic field 

Since the predecessor of the accretion disk is a highly- 
magnetic NS, the conservation of magnetic flux requires 
an inherited magnetic field in the disk. The successive 
magnetic flux from the remanent NS can be given by 



$ NS = 2tt 



/■r ou t,i 

>'l-in,l 



Birdr, 



(7) 



where B\ is the magnetic induction strength of the disk 
in the initial condition. The instabilities become effective 
when the magnetic pressure in the radial direction can 
support against gravity of the BH (e.g., Spruit & Uzden- 
sky 2005; Proga & Zhang 2006). The accretion process 
may be closed, if the magnetic induction strength satis- 
fies a critical value -B C rit,n, 



f>2 

^crit.n 
47T 



GM n £ n | 



(8) 



where M n and E n are the mass of the BH and surface 
density in the nth episode. The critical radius r cr ;t.n(_ 
Tout) can be estimated with 



B n rdr 



^crit,n?*i n)n ) 



(9) 



where B n is calculated by the conservation of magnetic 
flux. Therefore, the timescale is obtained by 



rp 1 n 

\Vn\ 



(10) 



where |i> n | is the absolute value of the average radial ve- 
locity from r cr it, n to 
We estimate the magnetospheric radius from 

r m « 6xl0 3 (v/v f f ) 2 / 3 (M/M Q s-^-^M/SM©)- 4 / 3 ^, 

(11) 

where v, Vff, M and r g (= 2GM/c 2 ) are the radial veloc- 
ity of the accretion disk, the free fall velocity, the mass 
accretion rate and Schwarzchild radius, respectively (e.g., 
Narayan et al. 2003; Proga & Zhang 2006). Thus the 
timescale of magnetic field dissipation can be estimated 
as ~ r m /vjf. We simplify the magnetic field as a uniform 
field with a magnetic flux 10 29 G cm 2 in our calculation 
(e.g., Proga & Zhang 2006). For Mi = 3 M Q , a»i = 0.9, 
Mi = 0.05 Mq s"\ ai = 0.01 and v/v ff ~ 10~ 2 - 10~ 3 , 
the timescale of the first episode is T\ ~ 2 s from Eq. 
(10), and the timescale of magnetic field dissipation is 
~ 0.1 — Is. If the remanent magnetic flux is reduced to 
5 x 10 28 G cm 2 and M 2 ~ Mi, a 2 ~ ai, a* 2 ~ a *i, the 
timescale of the second episode is T 2 ~ 27 s in case of 
M 2 = 0.01 Mq s _1 . Therefore, our model can potentially 
explain the EE of SGRBs. 

2.3. NDAF model 

We adopt the method present by Riffert & Herold 
(1995). This method is dedicated to numerically investi- 
gate the NDAF in vicinity of a rotating BH. The method 
defines general relativistic correction factors quoted as 



below, 
A = 



2GM 



+ 



5 = 1 - 



3GM 



/ a,GM \ 2 
V c 2 r / ' 



2a 



C = l-4a 



a.GMy 



c^r 



D = 



X c 



3xMc z 



4(: 



il/'' 



3aiM A 



G 



(xr) 2 



3xMc 2 



+ 2 (— L G ) 2 1 



(12) 
(13) 
(14) 

-cfa(15) 



where M, a* and r ms are the mass, dimensionless spin 
parameter and the radius of marginally stable orbit of the 
BH (e.g., Kato et al. 2008), respectively. The continuity 
equation remains valid, 



M 



-2nrl]v. 



(16) 



The hydrostatic equilibrium in the vertical direction 
leads to a corrected expression for the half thickness of 
the disk (Riffert & Herold 1995; Lei et al. 2009; Liu et 
al. 2010b), 

where c s = (p/ p) 1 ^ 2 is the isothermal sound speed, p and 
p are the total pressure and density of the disk, respec- 
tively. The viscous shear T r< ^ is also corrected as (Liu et 
al. 2010b) 

T« = -a Pj ^, (18) 

where a is a dimensionless constant that absorbs all the 
detailed microphysics of the viscous processes. The an- 
gular momentum equation can be simplified as 



M (GM\h/D\h 



\ r 3 ) 



(19) 



The total pressure includes the gas pressure from nu- 
cleons p gas , radiation pressure of photons p ra d, degener- 
acy pressure of electrons p c , and radiation pressure of 
neutrinos p v (see e.g., Liu et al. 2007), 



V = Pgas + Prad + Pe + Pv 

The energy equation is given by 

Qvis — Qadv "f" Qphoto "I" Qvi 



(20) 



(21) 



where <5 v i s , Qadv, Qphoto and Q v are the viscous heating 
rate, the advective cooling rate, the cooling rate due to 
photodisintegration of a-particles and the cooling due to 
the neutrino radiation, respectively (see e.g., Liu et al. 
2007). The heating rate Q v is is expressed as 



3GMM D 
87rr 3 B' 



(22) 



The radiation luminosity of the neutrinos released from 
the disk is obtained with the neutrino cooling rate Q„, 
i.e., 

(Tout 

L u = 4tt / Q v rdr, (23) 



/To, 

:„ = 4tt / 

J rin 
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We follow the approach by Ruffert et al. (1997), 
Popham et al. (1999), and Rosswog et al. (2003) to cal- 
culate the neutrino annihilation luminosity. The disk is 
modeled as a grid of cells in the equatorial plane. A cell 
k has its mean neutrino energy neutrino radiation 

luminosity l 1 ^., and distance to a space point above (or 
below) the disk d k - The angle at which neutrinos from 
cell k encounter antineutrinos from another cell k' at that 
point is denoted as d k k' ■ Then the neutrino annihilation 
luminosity at that point is given by the summation over 
all pairs of cells, 




jk lk f pk _i_ ^k! 

k u k k , u k' e v^Vi 

where A\ w 1.7 x 10~ 44 cm erg~ 2 s -1 and A 2 ~ 1.6 x 
10~ 56 cm erg" 2 s^ 1 (e.g., Popham et al. 1999). The 
total neutrino annihilation luminosity is integrated over 
the whole space outside the BH and disk, 

/>oo />oo 

L vV = Any2 / l u j7jdrdz. (25) 

However, as shown by Popham et al. (1999) and Liu 
et al. (2007), the neutrino annihilation would inject a 
highly beaming outflow around the inner part of the disk. 
For M = 3 Mq, M = 0.01 - 1 M© s" 1 and a = 0.01 ~ 
0.1, the opening angle of the ejection 9 is about 10° ~ 
20°. We conservatively assume that the opening angle 
of ejection is 10° and the efficiency of fireball is 77 = 0.1 
in our calculations. The observed isotropic luminosity of 
the nth step Li so , n can be estimated with 

iiso,n = T}L vii>n /{l - COS0). (26) 

3. NUMERICAL RESULTS 

Observationally, half of the Type I GRBs have EE 
component detection. We show the lightcurves of some 
SGRBs with EE component in Figure 2. The lightcurves 
are visually recognized as different emission episodes. We 
apply our model to these SGRBs. For simplicity, we 
consider only two emission episodes, i.e., the initial hard 
spikes and the EE component. We describe our numeri- 
cal method for two emission episodes as follows. 

There are seven unknown variables in our model, i.e., 
a*!, a» 2 , Mi, Mi, M 2 , Ml, and M 2 . Therefore, a group 
of seven equations are required in our calculations. Three 
of them are from the conservation of mass (energy), i.e., 

M 2 = Mi + MiT iein , n , (27) 

M 2 *=M*-M 1 Ti, (28) 

and 

M% = M 2 T 2 + SM* , (29) 

where 5M* is the mass of the residual disk post the sec- 
ond episode. Two other equations are from the conser- 
vation of angular momentum, i.e., 

J 2 = Ji + MiTil inA (30) 



and 

J 2 * = Jl - MiTil^i, (31) 

The values of J\ and J 2 are calculated with Eq. (4). 
Jl (J 2 *) is determined by Mi(M 2 ), a*i(a* 2 ), Mi(M 2 ), 
M{(M%), and a, which is calculated by Eq. (5). The 
rest two equations are related to L v „, which is a function 
of M, a*, M, and a in the NDAF model, 

L V u,i = f(Mi,a*i,Mi,a), (32) 

L vf >,2 = f(M 2 ,a* 2 ,M 2 ,a), (33) 

where the function / can be obtained from Eqs. (12)- 
(25). The average neutrino annihilation luminosity in 
each episode is calculated assuming that the accretion 
rate and spin parameter of the BH are a constant. The 
accretion timescale in each step depends on the initial 
mass of the disk and the switch of the magnetic field. The 
procedure of our calculations is descried as following. 

First, we assign the initial parameters of the BH and 
disk. We fix the initial mass and the spin parameter 
of the BH as 3Af Q and 0.9, respectively. The viscous 
parameter is assumed to be a = 0.01. The initial mass 
of the disk is adjustable in our calculation. 

Second, we take the observed average luminosity and 
timescale of the initial hard spikes as L v p t i and T\, then 
calculate M x with Eq. (32). The values of M 2 , M|, a* 2 , 
and M 2 are derived from Eqs. (27), (28), (30), and (33), 
respectively. 

Third, we calculate T 2 and L„p i2 with Eqs. (6) and (33) 
by taking a small value of 8M* , i.e., about ~ O.1M . 

We adjust the initial mass of the disk and find that in 
case of M* = 0.8Af© our results are roughly consistent 
with the observations for some typical SGRBs with EE, 
i.e., 050724, 060614, 061006, 061210, 070714B, 071227. 
Our results are shown in Figure 3. The first episodes of 
these GRBs are in the left-top circle of the figure. We 
correspondingly derive a region of the second episodes 
for the circle with the universal parameters mentioned 
above. We notice that besides GRB 050724 and 071227, 
the second episodes of the other four GRBs are in the 
region. Note that the EE component of the two GRBs 
are likely a late flare, which may have different physical 
origin as mentioned in §1. 

The mass of the disk is a very important factor for the 
EE. The decrease of M* would result in the significant 
decrease of the timescale and luminosity in the second 
episodes. For the case of Ml = 0.2Af Q , we find that the 
timescale of the second episode is shorter than 20 sec- 
ondf03 and the luminosity is roughly 1 order of magni- 
tude lower than that of the case Ml = 0.8M Q , as shown 
in Figure 3. Our results indicate that a bright EE com- 
ponent may be only detected for SGRBs with a massive 
disk. 

4. CONCLUSIONS AND DISCUSSION 

We have proposed that both the outward angular mo- 
mentum transfer and the switch of the magnetic barrier 
of a rotating BH-neutrino-cooled disk system may result 
in a long-lasting, impulse engine to produce several radi- 
ation episodes as observed in Type I GRBs. Based on the 

11 A low viscosity would increase the accretion time (e.g., Met- 
zger et al. 2008) . In our calculations, we take a\ = a.2 = 0.01 
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NDAF model, we have presented detailed calculations for 
some typical GRBs with detection of the EE component. 
Our numerical results in reasonable parameter sets well 
agree with the data. 

Observationally, about half of the Type I GRBs have 
EE component detection with Swift/BAT. This may be 
due to both the physical and instrumental effects. Our 
result suggests that a highly magnetic, massive disk is 
required to produce an intense, long-lasting EE compo- 
nent. Massive disk may from coalescence of a BH with 
a massive NS, i.e., M NS ~ 2 Mq (Morrison et al. 2004; 
Dai et al. 2006; Demorest et al. 2010; Li et al. 2012). 
Kluzniak & Lee (1998) showed that the mass of the disk 
from merger of a BH-NS binary may be larger than 0.5 
Mq (see also Janka et al. 1999). For typical GRBs with 
EE detection by Swift/BAT, our model suggests that the 
disk should be ~ 0.8 Mq. In case of the mass being lower 
than ~ 0.2 Mq, the EE may cannot be detected with 
Swift/BAT and its timescale is shorter than 20 seconds. 
Both analytical and simulated investigations for coales- 
cence of compact objects have been extensively studied 
(e.g., Ruffcrt et al. 1997; Kluzniak & Lee 1998; Lee & 
Kluzniak 1999; Ruffert 1999; Lee & Ramirez-Ruiz 2002; 
Lee et al. 2005; Proga & Zhang 2006; Ruffert & Janka 
2010). The suggested mass of the disk by these authors 
is much lower than that predicted by our model. Con- 
sidering that the mass of a NS is 1 — 1.4 Mq, Ruffert & 
Janka (1997) pointed out that the coalescence of the NS- 
NS binary might result in a BH ~ 2.5 Mq surrounding 
by a disk with mass ~ 0.1 — 0.2 Mq. The key ingre- 
dients in our model are a massive disk and the angular 
momentum transfer, which may significantly prolong the 
lifetime of the GRB central engine for the EE compo- 
nent. Therefore, the EE component may be a probe for 
merger process and formation of massive disk in the cen- 
tral engine of SGRBs. 

The detection or not of the EE component might be 
also due to instrumental selection effect. Our model sug- 
gests that a disk with M* < 0.2M© would not produce an 
EE component that can be detectable with Swift/BAT- 
like instruments. Assuming a constant radiation effi- 
ciency, the observed luminosity is proportional to the 
neutrino annihilation luminosity. Since the peak energy 
of the vf v spectrum is tightly correlated with the radia- 
tion luminosity (Liang et al. 2004; Yonetoku et al. 2004), 
the EE component of some SGRBs may be out of the in- 
strument bands and only the initial hard-short episode 
can be detectable as usually seen in typical SGRBs. 
With the results for M{ = 0.8A/© and' M{ = O.2M 
one can observed a luminosity-duration correlation for 
the EE component. This seems to be at odds with ob- 
servations that show the opposite trend: longer events 
have a lower luminosity. Note that the observations for 
the EE component greatly suffer from the selection ef- 
fect of instrumental sensitivity. We check the observed 
luminosity-duration correlation for both the initial hard 



Bardeen, J. M. 1970, Nature, 226, 64 

12 The jet opening angle of SGRBs may be even much smaller 
than 10°, such as ~ 0.3° for GRB 090510 (Corsi et al. 2010; He et 
al. 2011; Fan & Wei 2011) 



spikes and the EE components for the GRBs with EE 
detection reported in Zhang et al. (2009), but no statisti- 
cally accepted correlation is found. The relation between 
the luminosity and duration of our model is parameter- 
dependent. This calls for a detailed analysis based on 
Monte Carlo simulations in order to elaborate this rela- 
tion and its observational biases. 

Some caveats for our model should be discussed. The 
neutrino annihilation luminosity is correlated positively 
with the accretion rate. In order to ensure the mass of 
NS less than about 2 Mq, a moderate accretion rate is 
required in our model. The derived accretion rates for 
both the initial and late episodes from our model are less 
than 0.05 Mq s _1 , which is lower than the typical ones, 
i.e., 0.1 ~ IMq s . Note that the typical values are 
for a spatially uniform neutrino annihilation luminosity. 
In our calculations, the spatial distribution of the neu- 
trino annihilation luminosity is collimated (e.g., Popham 
et al. 1999; Liu et al. 2007). The jet opening angle is 
conservatively assumed to be 10° in our calculation^. 
Therefore, the isotropic-equivalent accretion rate would 
be larger than 1^2 order of magnitude than our re- 
sults, which is consistent with the typical values used by 
previous authors. 

We have considered the switch of magnetic field for 
the accretion processes, but the effects to confine the 
jet opening angle by the BZ mechanism (Blandford & 
Znajek 1977) and the enhance of the neutrino annihi- 
lation are ignored. These effects may also greatly im- 
pact on the jet luminosity. Lei et al. (2009) investi- 
gated the magnetic coupling between the BH and disk. 
They found that the luminosity of neutrino annihilation 
is much larger than the luminosity of NDAF without 
magnetic field. Recently, Barkov & Pozanenko (2011) 
proposed a two jet model which describes both main 
and EE components by different off-axis position of ob- 
server. Their model involves a short duration jet powered 
by heating due to neutrino annihilation and a long-lived 
BZ jet with significantly narrow opening angle. The BZ 
mechanism can replace the neutrino annihilation to pro- 
duce the main emission and EE in our model. Since the 
lack of information of the intensity and distribution of 
the magnetic field, we restrict to discuss the neutrino 
annihilation as the main power source in our model. 

ACKNOWLEDGMENTS 

We thank Wlodzimierz Kluzniak, Bing Zhang, Feng 
Yuan, Ye-Fei Yuan, Yi-Zhong Fan, Li Xue, Xiao-Hong 
Zhao, and Da-Bin Lin for beneficial discussions and 
comments. This work is partially supported by the 
National Basic Research Program ("973" Program) of 
China (Grant 2009CB824800), the National Natural Sci- 
ence Foundation of China (Grants 10873005, 10873009, 
11003004, 11025313, 11033002, 11073015, 11103015, 
11222328 and 11233006), and Guangxi Science Foun- 
dation (2010GXNSFC013011 and special support with 
Contract No. 2011-135). 



Barkov, M. V., & Pozanenko, A. S. 2011, MNRAS, 417, 2161 
Barthelmy, S. D., Cannizzo, J. K., Gehrels, N., et al. 2005, ApJ, 
635, L133 

Belczynski, K., Holz, D. E., Fryer, C. L., et al. 2010, ApJ, 708, 117 
Blandford, R, D. & Znajek, R. L. 1977, MNRAS, 179, 433 



6 



Liu et al. 



Burrows, D. N., Grupe, D., Capalbi, M., et al. 2006, ApJ, 653, 468 
Cao, X. 2011, ApJ, 737, 94 

Chen, W. X., & Beloborodov, A. M. 2007, ApJ, 657, 383 
Connaughton, V. 2002, ApJ, 567, 1028 
Corsi, A., Guetta, D., & Piro, L. 2010, ApJ, 720, 1008 
Dai, Z. G, & Lu, T. 1998, A&A, 333, L87 

Dai, Z. G, Wang, X. Y., Wu, X. F., & Zhang, B. 2006, Science, 
311, 1127 

Delia Valle, M., Chincarini, G., Panagia, N., et al. 2006, Nature, 
444, 1050 

Demorest, P. B.,Pennucci, T., Ransom, S. M., Roberts, M. S. E., 

& Hessels, J. W. T. 2010, Nature, 467, 1081 
Di Matteo, T., Perna, R., & Narayan, R. 2002, ApJ, 579, 706 
Donaghy, T . Q., Lamb, D. Q., Sakamoto, T., et al. 2006, 

|arXiv:astr o-ph/0605570 
Eichler, D., Livio, M., Piran, T., & Schramm, D. N. 1989, Nature, 

340, 126 

Fan, Y.-Z., & Wei, D.-M. 2005, MNRAS, 364, 42 
Fan, Y.-Z., & Wei, D.-M. 2011, ApJ, 739, 47 
Fynbo, J. P. U., Watson, D., Thone, C. C, et al. 2006, Nature, 
444, 1047 

Gal- Yam, A., Fox, D. B., Price, P. A., et al. 2006, Nature, 444, 
1053 

Gehrels, N., Norris, J. P., Barthelmy, S. D., et al. 2006, Nature, 
444, 1044 

Goldstein, A., Preece, R. D., & Briggs, M. S. 2010, ApJ, 721, 1329 

Gu, W.-M., Liu, T., & Lu, J.-F. 2006, ApJ, 643, L87 

He, H.-N., Wu, X.-F., Toma, K., Wang, X.-Y., & Meszaros, P. 

2011, ApJ, 733, 22 
Janka, H.-T., Eberl, T., Ruffert, M., & Fryer, C. L. 1999, ApJ, 

527, L39 

Janiuk, A., Moderski, R., & Proga, D. 2008, ApJ, 687, 433 
Kato, S., Fukue, J., & Mineshige, S. 2008, Black-Hole Accretion 

Disks: Towards a New Paradigm (Kyoto: Kyoto Univ. Press) 
King, A., O'Brien, P. T., Goad, M. R., et al. 2005, ApJ, 630, L113 
Kluzniak, W., & Lee, W. H. 1998, ApJ, 494, L53 
Kohri, K., & Mineshige, S. 2002, ApJ, 577, 311 
Kohri, K., Narayan, R., & Piran, T. 2005, ApJ, 629, 341 
Kouveliotou, C, Meegan, C. A., Fishman, G. J., et al. 1993, ApJ, 

413, L101 

Lazzati, D., Blackwell, C. H., Morsony, B. J., & Begelman, M. C. 

2011, MNRAS, 411, L16 
Lazzati, D., Perna, R., & Begelman, M. C. 2008, MNRAS, 388, 

L15 

Lazzati, D., Ramirez-Ruiz, E., & Ghisellini, G. 2001, A&A, 379, 
L39 

Lee, W. H., & Kluzniak, W. 1999, ApJ, 526, 178 
Lee, W. H., & Ramirez-Ruiz, E. 2002, ApJ, 577, 893 
Lee, W. H., & Ramirez-Ruiz, E. 2007, New Journal of Physics, 9, 
17 

Lee, W. H., Ramirez-Ruiz, E., & Lopez-Camara, D. 2009, ApJ, 
699, L93 

Lee, W. H., Ramirez-Ruiz, E., & Page, D. 2005, ApJ, 632, 421 
Lei, W. H., Wang, D. X., Zhang, L., et al. 2009, ApJ, 700, 1970 
Levesque, E. M., Bloom, J. S., Butler, N. R., et al. 2010, 

MNRAS, 401, 963 
Li, A., Huang, F., & Xu, R.-X. 2012, Astroparticle Physics, 37, 70 
Liang, E. W., Dai, Z. G., & Wu, X. F. 2004, ApJ, 606, L29 
Lin, L., Liang, E.-W., Zhang, B.-B., & Zhang, S. N. 2008, AIPC, 

1065, 39 



Lin, L., Liang, E.-W., & Zhang, S. N. 2010, Science in China G: 

Physics and Astronomy, 53, 64 
Liu, T., Gu, W.-M., Dai, Z.-G, & Lu, J.-F. 2010a, ApJ, 709, 851 
Liu, T., Gu, W.-M., Xue, L., & Lu J.-F. 2007, ApJ, 661, 1025 
Liu, T., Gu, W.-M., Xue, L.. & Lu J.-F. 2012, Ap&SS, 337, 711 
Liu, T., Gu, W.-M., Xue, L., Weng, S.-S., & Lu J.-F. 2008, ApJ, 

676, 545 

Liu, T., Liang, E.-W., Gu, W.-M., et al. 2010, A&A, 516, A16 
Lii, H.-J., Liang, E.-W., Zhang, B.-B., & Zhang, B. 2010, ApJ, 
725, 1965 

Lyons, N., O'Brien, P. T., Zhang, B., et al. 2010, MNRAS, 402, 
705 

Mctzger, B. D., Arcones, A., Quataert, E., & Martmez-Pinedo, G. 

2010, MNRAS, 402, 2771 
Metzger, B. D., Quataert, E., & Thompson, T. A. 2008, MNRAS, 

385, 1455 

Morrison, I. A., Baumgarte, T. W., & Shapiro, S. L. 2004, ApJ, 
610, 941 

Nakar, E. 2007, Phys. Rep., 442, 166 

Narayan, R. , Igumenshchev, I. V., & Abramowicz, M. A. 2003, 
PASJ, 55, L69 

Narayan, R., Paczynski, B., & Piran, T. 1992, ApJ, 395, L83 
Narayan, R., Piran, T., & Kumar, P. 2001, ApJ, 557, 949 
Norris, J. P., Gehrels, N., & Scargle, J. D. 2010, ApJ, 717, 411 
Paczynski, B. 1991, Acta Astronomica, 41, 257 
Paczynski, B. 1998, ApJ, 494, L45 
Perley, D. A., et al. 2009, ApJ, 696, 1871 

Perna, R., Armitage, P. J., & Zhang, B. 2006, ApJ, 636, L29 

Piran, T. 2004, Rev. Mod. Phys., 76, 1143 

Popham, R., Woosley, S. E., & Fryer, C. 1999, ApJ, 518, 356 

Proga, D., & Zhang, B. 2006, MNRAS, 370, L61 

Qin, Y., Liang, E.-W., Yi, S.-X., et al. 2012. larXiv: 1205.11881 

Riffcrt, H., & Herold, H. 1995, ApJ, 450, 508 

Rosswog, S. 2007, MNRAS, 376, L48 

Rosswog, S., Ramirez-Ruiz, E., & Davies, M. B. 2003, MNRAS, 

345, 1077 
Ruffert, M. 1999, A&A, 346, 861 

Ruffert, M., & Janka, H.-T. 1997, Reviews in Modern Astronomy, 
10, 201 

Ruffert, M., k Janka, H.-T. 2010, A&A, 514, 66 

Ruffert, M., Janka, H.-T., Takahashi, K., & Schaefer, G. 1997, 

A&A, 319, 122 
Spruit, H. C, & Uzdensky, D. A. 2005, ApJ, 629, 960 
Sun, M.-Y., Liu, T., Gu, W.-M., & Lu, J.-F. 2012, ApJ, 752, 31 
Thorne, K. S. 1974, ApJ, 191, 507 

Virgili, F. J., Zhang, B., O'Brien, P., & Troja, E. 2011, ApJ, 727, 
109 

Wang, D. X., Xiao, K., & Lei, W. H. 2002, MNRAS, 335, 655 
Woosley, S. E. 1993, ApJ, 405, 273 
Woosley, S. E., & Bloom, J. S. 2006, ARA&A, 44, 507 
Xin, L.-P., Liang, E.-W., Wei, J.-Y., et al. 2011, MNRAS, 410, 27 
Yonetoku, D., Murakami, T., Nakamura, T., et al. 2004, ApJ, 
609, 935 

Yuan, F., & Zhang, B. 2012, ApJ, 757, 56 

Zhang, B. 2006, Nature, 444, 1010 

Zhang, B., & Meszaros, P. 2002, ApJ, 566, 712 

Zhang, B., & Meszaros, P. 2004, Int. J. Mod. Phys. A, 19, 2385 

Zhang, B., Zhang, B.-B., Liang, E.-W., et al. 2007, ApJ, 655, L25 

Zhang, B., Zhang, B.-B., Virgili, F. J., et al. 2009, ApJ, 703, 1696 



Radial Angular Momentum Transfer and Magnetic Barrier for SGRB Central Engine Activity 



7 




Fig. 1. — Schematic illustration of our model and corresponding cartoon lightcurvc: (a) the initial state of the central engine — the filled 
circle stands for the BH and the fuscous trapeziform region for the disk with magnetic field (curves); (fe n ) angular momentum transfer 
process and outward flow (light gray region); (c n ) magnetic barrier in vicinity of the BH, where n is for the nth emission episode. 
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Fig. 2. — BAT lightcurves of six SGRBs with detection of the EE component. 
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Fig. 3. — Our numerical results for the six SGRBs shown in Figure 2. The filled symbols and open symbols figure the main and the 
corresponding EE components of the SGRBs. The first episodes of these GRBs are in the left-top circle of the figure. We correspondingly 
derive a region of the second episodes with the universal parameters for M-j* = 0.8Mq. Besides GRB 050724 and 071227, the second 
episodes of the other four GRBs are in the region. The EE component of the two GRBs are likely a late flare, which may have different 
physical origin as mentioned in §1. The corresponding region for M£ = O.2M0 is also shown. It indicates that the luminosity of the second 
episode is roughly 1 order of magnitude lower than that for the case of Mj* = 0.8Mq and the timescale is shorter than 20 seconds. The 
switch timescale of the magnetic barrier is 0.1-1 seconds. This timescale is much smaller than the observed timescale of the EE component. 
We thus ignore it here. 



